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We explore analytically the nature of the transition to the Fulde-Ferrel- Larkin-Ovchinnikov 
superfluid phases in the vicinity of the tricritical point, where these phases begin to appear. We 
make use of an expansion of the free energy up to an overall sixth order, both in order parameter 
amplitude and in wavevector. We first explore the minimization of this free energy within a subspace, 
made of arbitrary superpositions of plane waves with wavevectors of different orientations but same 
modulus. We show that the standard second order FFLO phase transition is unstable and that 
a first order transition occurs at higher temperature. Within this subspace we prove that it is 
favorable to have a real order parameter and that, among these states, those with the smallest 
number of plane waves are prefered. This leads to an order parameter with a cos(qo.r) dependence, 
in agreement with preceding work. Finally we show that the order parameter at the transition is 
only very slightly modified by higher harmonics contributions when the constraint of working within 
the above subspace is released. 

PACS numbers : 74.20.Fg, 74.60. Ec 



I. INTRODUCTION 

Although they have been proposed a long time ago, Fulde-Ferrel- Larkin-Ovchinnikov (FFLO) phases are still 
the subject of a continuing interest. Indeed the existence of these phases is a fairly remarkable phenomenon since 
they correspond to a spontaneous symmetry breaking of the standard BCS superfluid phase in the presence of an 
effective field, inducing a difference in chemical potential between the two populations involved in the formation of 
Cooper pairs. This symmetry breaking leads to an inhomogeneous superfluid with a space dependent order parameter, 
while the applied field is perfectly homogeneous. This situation is analogous to the appearance of vorticity in type II 
superconductors, but in this latter case the effect is due to the coupling of the field to particle currents while in FFLO 
phases only the coupling to the spin of the pairing fermions is involved. In standard superconductors the coupling to 
the orbital degrees of freedom is much stronger than the coupling to the spins. Hence the upper critical field is due to 
the orbital coupling and the FFLO phases can not be observed, since they should appear at much higher field. However 
in heavy fermions superconductors the strength of these two couplings is comparable, which could make possible the 
observation of FFLO phases. Nevertheless their sensitivity to impurities could be a major problem. Another possible 
direction to eliminate the orbital coupling is to consider lower dimensional superconductors, in a geometry where the 
currents would have to flow in an actually prohibited direction. Organic compounds or cuprate superconductors are 
interesting systems in this respect. And indeed very recently the FFLO state has been claimed to be observed in 
a quasi-two-dimensional organic compound |^ . On the other hand earlier possible observations in heavy fermion 
compounds have not been undisputed. We note in particular that the analysis of experimental results relies very 
often heavily on the theoretical results, but we will see that the situation is not completely satisfactory in this respect. 

Another class of physical systems where FFLO phases could be observed is coming up quite recently. These are 
^ , the ultracold fermionic gases. As it is well known remarkably low temperatures have been obtained on bosonic 
k>( ' gases, leading in particular to the observation of Bose-Einstein condensation in alkali ultracold gases. More recently 
r> ' fermionic gases have been cooled down in the degenerate regime ||^-|^ and reaching a BCS superfluid transition in 
^ these systems seems a reasonable possibility [^,^ . However in the systems considered for observing this transition 
there is, in contrast to electronic spin relaxation in superconductors, no fast relaxation mechanism to equalize the 
populations of the two fermion species involved in the formation of Cooper pairs. Hence one should have no limitation 
to the effective field in these systems since the number of atoms in the populations can be in principle obtained at 
will. So the difference in atomic populations looks as a very promising control parameter. On the other hand if this 
parameter is not fully controlled this might very well be a major difficulty in reaching the BCS transition in these 
systems . 

In contrast with this raising experimental interest there are still theoretical problems with the precise nature of 
the possible phases. Specifically the basic FFLO instability corrresponds to have pairs formed with a total nonzero 
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momentum qo instead of forming pairs (k, — k) with zero total momentum as in the standard BCS phase. This 
gives rise to a spatial dependence exp(iqo.r) for the order parameter, which leaves a degeneracy with respect to the 
orientation of qo. This has been investigated by |2| Larkin and Ovchinnikov (LO) who looked how it is lifted right 
below the critical field. In this case, when considering the spatial dependence of the order parameter A(r), one 
can restrict the investigation to the subspace generated by linear combination of the plane waves exp(zqo.r) with 
all possible directions qo. At T = 0, LO looked for periodic structures and found that the energetically favored 
result is a second order transition to a one-dimensional 'planar' texture A(r) ^ cos(qo.r). However they left open in 
their paper the possibility of a first order transition. Actually when considering the three-dimensional 'cubic' texture 
A(r) ^ cos(goa;) + cos{qoy) + cos{qoz) they found that it is energetically unfavorable compared to the normal state. 
Nevertheless they obtained from the gap equation that a nonzero solution for this order parameter exists above the 
FFLO transition line. In terms of the expansion of the free energy in powers of the order parameter, schematically 
f2 = (22^^ -|- a4A'*, this situation corresponds to a positive coefficient a2 for the second order term and a negative 
coefficient for the fourth order one, just the opposite of the standard Landau- Ginzburg expansion below the 
transition. The LO evaluation of the related free energy corresponds actually to the maximum Q!2/4|a4| of this free 
energy. Beyond this maximum the free energy decreases, and it would go to —oo if one would consider only the second 
and fourth order terms. Naturally one has to include the effect of all higher order terms in order to find the value 
of the free energy for large values of the order parameter. However, at the FFLO transition line and slightly above 
it, the free energy becomes negative for values of the order parameter (specifically for A'^ = a2/|a4|) where only the 
second and fourth order terms have to be kept in the expansion. This shows definitely that the transition toward the 
FFLO phase occurs above the standard second order FFLO transition line and is actually first order, because one 
can display in this range a solution which has a lower free energy than the normal state. On the other hand in order 
to obtain consistently the order parameter which gives the lowest free energy one has naturally to take into account 
higher order terms. In particular it is by no means obvious that the cubic phase is the stable one. Moreover when 
higher order terms are considered there are no reason anymore to restrict the search to the subspace generated by the 
plane waves exp(iqo.r). 

In order to explore more fully this difficult problem with easier conditions, it is better to be able to proceed to some 
kind of expansion. This can be done if, instead of working at T = 0, one explores the vicinity of the tricritical point 
(TCP), where the FFLO transition line starts. It is located at Ttcp/Tco = 0.561 where Tco is the critical temperature 
for p, = , with 2/1 = /X| — /i| being half the chemical potential difference between the two fermionic populations 
forming pairs. The corresponding effective field is fttcp/Tco = 1.073. At this point, in the free energy expansion, both 
the coefficient of the second ordre term a2 and the coefficient of the fourth ordre term vanish. Naturally a2 is 
zero just because the TCP is on the standard second order phase transition line. On the other hand following this 
transition line one has Q!4 > for T > Ttcp and 0:4 < for T < Ttcp. The change of sign of at Ttcp is the origin 
of the FFLO instability. Clearly, by continuity, the equilibrium order parameter will be small in the vicinity of the 
TCP since we know that it is zero just above the TCP on the second order line and it is small right below it in the 
superfluid phase. Therefore a power expansion of the free energy will be enough to find it. On the other hand, since 
the second and fourth order terms are zero at the TCP, we have clearly to expand at least up to sixth order, but 
this will prove to be enough because the corresponding coefficient is positive and not small. Similarly the optimum 
wavevector qo corresponding to the FFLO phase will be small in the vicinity of the TCP since it is zero just above it 
on the second order line. This allows to proceed to a gradient expansion of the free energy. Since only even powers 
of the wavevector can enter and we look for a minimum as a function of this wavevector, we have to expand at least 
up to fourth order in gradient, but this will prove again to be enough. 

This line of thought has actually already been followed by Houzet et al. |l^Jl^ , who have performed this expansion 
for the free energy and explored the result numerically. They have found that the energetically favored phase is the 
one-dimensional planar order parameter found by LO at T = 0, but that the transition is actually slightly first order, 
instead of second order as found by LO at T = 0. Our purpose in the present paper is rather to proceed to an 
analytical study of this problem. Indeed there are infinitely many possible order parameters in competition. And our 
aim, in considering the vicinity of the TCP, is to find the important ingredients which are responsible for the selection 
of the actual stable state and obtain a better physical understanding, having in particular in mind the generalization 
to more complicated situations. Hence our paper is complementary to their work. In particular we obtain a first order 
transition to the one- dimensional planar order parameter, but we will be able to analyze the reasons which favor 
this phase. The transition to this planar order parameter has been actually explored numerically down to T = by 
Matsuo et al. . They have used quasiclassical equations and found that the transition keeps first order down to 
low temperature, but eventually goes to second order in agreement with LO. On the other hand since we know that 
at T = the cubic phase is more stable than the planar one, the question of the stablest phase at low temperature is 
still unsolved. Finally this first order transition to the planar phase in the three-dimensional case is in contrast with 
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the results of Burkhardt and Rainer |14| who found it to be second order in a two-dimensional space. 



In the following section, for completeness and to set up our notations, we rederive the expression pq ] of the free 
energy. After considering in section III some simple situations, we explore in details the minimization of this free 
energy. This is done in section IV by restricting our search to the LO subspace for the order parameter, which is made 
of arbitrary superpositions of plane waves. In section V we show that our results are only slightly modified when 
we release this restriction. Throughout the paper we restrict ourselves to the simplest BCS scheme, namely we will 
consider the free energy corresponding to a weak coupling isotropic Fermi system, ignoring in particular any Fermi 
liquid effect. Moreover we concentrate on the three-dimensional case which leads to a first order transition, and make 
only occasionally comparison with the two-dimensional situation where the transition is second order. 



II. THE FREE ENERGY 



There are various ways to obtain the explicit expression for the free energy we need in the vicinity of the TCP 
[^,0. In practice it is convenient to use the fact jl^ that, by varying the free energy with respect to A(r) one finds 
the gap equation, which is also easily obtained from Gorkovs equations, as done for example by LO |^ . The integral 
form of these equations is , with standard notations: 

Gt(r,r') = G?(r-r')- / driG?(r - ri)A(ri)F+(ri, r') (1) 



^^+(r,r') = y driG;(r-ri)A*(ri)GT(ri,r') (2) 

with, for the Fourier transforms of the free fermions thermal propagators, G^(k) — [iojn — + A*)"^ and Gj'(k) = 
{—iujn — '^fe — m)^^ where is the kinetic energy measured from the Fermi surface for /i = and a;„ = ■nT{2n + 1) 
are Matsubara frequencies. The order parameter is given by the self- consistency relation : 

/^*{v) = VTY,F+{v,v) (3) 

n 

We expand Eq. (1-2) up to fifth order in A(r). We introduce the Fourier transform Aq = / (irA(r) exp(— iq.r) of the 
order parameter. As explained in the introduction we proceed also to an expansion in the wavevector q of the order 
parameter since we know that its relevant values will be small in the vicinity of the TCP. More precisely we will see 
that, in order to obtain a coherent expansion, it is enough to go only up to fifth order terms in overall power of A 
and q. This means that, in the gap equation, we have to expand the first order term in A only to fourth order in q. 
Similarly the third order term in A has to be expanded only to second order in q and the fifth order term in A can 
be calculated to zcroth order in q. For example in order to find the third order term, we have to expand up to second 
order in wavevectors : 

^ GO(k)G°(k + qi)GO(k -t- q - q3)G°(k + q)A;^ Aq, A^^ (4) 
k 

where we have used qi -I- qa = q + q2 and we have omitted the unnecessary spin index. In the expansion appear the 
following numerical coefficients: 

«o(/^,r) = ^-2^r3?E„.o^] 



a2WT) = -M22^r5RE„.o4] (5) 

a4(Ai/T) = -M*2^r5RE„=o4] 

where Qn = — and the summation for oq has to be cut-off in the standard BCS way. The simple second order 
transition line to a standard BCS superfiuid with space independent order parameter is given by ao[fi, T) = 0. Below 
the TCP it corresponds to a spinodal transition line, at which the normal state becomes absolutely unstable against 
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a transition toward a space independent order parameter. The domain ao(/i,T) > corresponds to the region of the 
(/I, T) phase diagram above this hne, and it is the domain where we will look for other transitions. In practice we can 
see aoifi, T) as a measure of the distance from the spinodal line in the (/i, T) plane. Explicitely if we define Tsp{ft/T) 
the spinodal temperature as a function of the ratio p,/T, we have ao{fi,T) = ln[T/Tsp{fL/T)]. We will not need to 
explicit further this distance. As indicated in the introduction we have by definition a2(/x/T) = at the TCP and it 
is small in the vicinity of this point. For (fi/T) > (/x/T)tcp = 1.91, we have 02 > and a2 < for {p./T) < {fi/T)tcp- 
Finally a4,{p,/T) = 0.114 at the TCP (while it is negative near Ji = Q and goes to -0.25 when T — » 0). With these 
notations the gap equation in the vicinity of the TCP reads : 

Aq[ao - \a2Q^ + \aiQ'\ - ^ A^,A*^^A^,[^a2 - ^a^iQ"" + 2Ql + 2Ql - Q.(Qi + Q3) + 3Q1.Q3)] 

qi 

+^«4^Aq,A;^Aq3A*^Aq,=0 (6) 
qi 

where we have used the dimensionless wavevector Q = qvp/^p and expressed Aq in units of p,. Also the momentum 

conservation is assumed in the summations, that is qi+qs = q+q2 in the third order term and qi+qa+qB = q+q2+q4 
in the fifth order one. The above expression can be checked against the case of the simple Fulde-Ferrell state 
A(r) = exp(iqo.r) where a single wavevector enters. 

Now the above gap equation is obtained by minimizing the following free energy difference J7 between the superfluid 
and the normal state : 

n = Y, |AqP[ao - ia2g2 + la^Q^] - 1 ^ Aq, A*^Aq3 A* jia2 - ^a^iQl + Ql + Ql + Ql + Q1.Q3 + Q2-Q4)] 
q qi 

+ia4^Aq,A;^Aq3A;^Aq,A;^ (7) 
qi 

where we have the momcntiim conservation qi + qa = q2 + q4 in the fourth order term while qi + qs + qs = q2 + q4 + qe 
holds in the sixth order one. We have used symmetry and momentum conservation to present the fourth order term 
in a symmetrical way. This expression Eq.(7) is just the free energy we were looking for. We have considered here 
the 3-D case. For a two- dimensional system the angular averages found in the calculation are different. The result 
is simply obtained from the above one by multiplying the terms by 3/2 and the Q** terms by 15/8. 

III. SIMPLE CASES 

Let us first consider some simple situations. If we consider an homogeneous order parameter, that is q = 0, we 
have merely O = aoA^ — a2A'*/4 + qaA^/S . If we want to have this free energy negative for ao > we need to 
have a2 > 0, that is to be at temperature below the TCP. In this case f2 > when ao > 02/804 and we reach a 
first order transition for ao = 02/8^4, with a non zero order parameter A^ = 02/04. This is the standard first order 
Pauli limiting transition. We consider next the possibility of a second order transition. In this case only the second 
order term in Eq.(7) is relevant. The location of this transition is given by oq = ^a2Q^ — ^04(5"*. Below the TCP, 
where 02 > 0, we can find oq > for non zero wavevector Q, that is we will find an FFLO phase. Precisely the 
optimal wavevector is Qq = §02/04 and the corresponding maximal oo is oq = ^02/04. We see that this value is 
larger than the one we just found for the standard Pauli limiting transition. Thus as expected the FFLO transition 
happens first and overtakes the first order transition. Finally it is natural and interesting to try to generalize the two 
above situations and consider the possibility of a first order transition for an order parameter with a single wavector 
component Aq. With the shorthand Aq = A the free energy writes : 

n = [ao- ^02Q^ + ^04Q^]A2 - l[o2 - 2o4Q^]A4 + ^04A6 (8) 

Minimizing first with respect to we obtain for the extremum the condition = §02/04 — |A^, which implies 
that must have A^ < I02/04 otherwise we are back to the homogeneous situation and the Pauli limiting transition. 
Inserting this value for Q"^ in Eq.(8) we find O = [oq — ^a2/a4]A^ + g02A^ — ^a4A^. We are naturally interested 
in finding a transition higher than the standard FFLO. This means we are looking for ao > ^02/04, so the first term 
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in is positive. But one sees that the sum of the last two terms is also positive for < |a2/a4. Therefore we have 
not been able to improve the standard FFLO solution. However we have clearly not done our best in this direction. 

Before trying to improve in this way, it is convenient to simplify our expression for the free energy by taking 
reduced units for the order parameter and the wavector, which come out naturally from our above discussion. We 
set A = (a2/a4)^/^A, Q = (a2/a4)^/^q , ao = Aoa\/ai and f2 = {a\/a\)F. This leads to rewrite Eq.(7) for the free 
energy as : 

F = Y, |Aq|^[Ao - \ct + \t] - E Aq3 A;^ [i - + ql + ql+ql + q^.qg + q^.q^)] 

q q* 

+^ENiA;.N3A;,Aq,A*, (9) 

q; 

It is clear from this rescaling transformation and the resulting expression that, in the vicinity of the TCP (where 02 
is small), it is unnecessary to go beyond our sixth order expansion in A and q. It is also of interest to rewrite this 
free energy as a functional of A(r) by Fourier transform. This gives, after by parts integrations : 

F = j dv[Ao\K\^ - 1|VA|2 + i|v2A|2] - j cir[l|A|4 - ^^[2{V\K\^f + 3(VA2)(VA*2)]] + j dv\K\^ (10) 



IV. THE LO SUBSPACE 

As emphasized by LO all the states corresponding to the same wavector go but with different orientation for qo 
are degenerate right on the FFLO transition line. This degeneracy is lifted, at least partially, when one goes into the 
superfiuid phase because of the coupling between the various plane waves produced by the nonlinear terms in the free 
energy. When one investigates the states selected in this process, one has to consider the subspace: 

A(r) = EAq„exp(iqo.r) (11) 

of all the order parameters generated by tlicjse plane waves. We call this the LO subspace. In this section we will 
restrict to this subspace our search for the state appearing at the transition : we will look for the minimum of the free 
energy within this LO subspace. Actually LO looked for a lattice as a solution and restricted themselves to this kind 
of order parameter. However there is physically no basic reason to enforce this type of restriction. One could look for 
incommensurate structures or quasicrystal-like solutions. Even if these are not the lowest energy solution, they might 
be of interest as local solutions corresponding physically to defects. Therefore we have not set a periodicity condition 
on the solutions we have considered. Nevertheless let us indicate at once that the energetically favored solutions we 
have found within the LO subspace are actually periodic. Although considering only the LO subspace is an important 
restriction, it does not make the problem easy at all, although we solve it completely below. 

Let us first show that, within this subspace, the standard FFLO transition line is not stable, and the transition is 
actually first order. With our reduced units the LO subspace corresponds to qg = |. This minimizes the coefficient of 
the second order term in Eq.(9). The FFLO transition line is then given by ^0 = which makes the second order 
term zero within the LO subspace. Let us then look at the fourth order term in Eq. (9) . This amounts to calculate the 
functional derivative of the free energy on the FFLO transition line. In a standard second order phase transition it 
should always be positive, forcing the order parameter to be zero on the transition line. However in the present case 
it is not obvious that this is systematically so because of the interplay of the wavevectors in this term. Specifically 
we introduce a parameter /? which describes this effect for any fixed order parameter A(r). It is defined by : 

2/3^2 ^Aq,A;^Aq3A;^ = g2^(qi.q3 +q2.q4)Aq,A*^Aq3A;^ = - / A2(VA*)2 +C.C. (12) 
qi qi ^ 

For the simple case of the FF solution we have merely /3 = 1. However this is the highest possible value and we can 
think of decreasing it, or even making it negative, by a proper choice of the order parameter in the LO space, although 
naturally j3 > —1. Then the fourth order term in Eq.(9) is given by [i(/3 + 2)^Q — 3] / |A|'' and has thus the same sign 
as 5/3 + 1 when we take = |- Hence for any order parameter with (3 < —\ the fourth order term is negative. We 
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will find actually many such states. Now when we arc on the FFLO transition line (the second order term is zero) and 
the fourth order term is negative, we decrease the free energy and make it negative just by taking the order parameter 
to be small and nonzero, which shows that the standard FFLO transition line is unstable. Naturally the larger the 
order parameter, the lower the free energy, but we have to stay in the range where the second and fourth order terms 
are the only ones important in our expansion, which means that the sixth order term is negligible. Then since we 
have a negative free energy, we can raise it back to zero by increasing beyond its value ^ on the FFLO line, 
which means we go beyond this line in the T) phase diagram. In this way the second order term becomes positive. 
Eventually we will be limited by the growth of the sixth order term. The transition we have found corresponds to 
positive second and sixth order terms and a negative fourth order one. The free energy becomes negative for a nonzero 
value of the order parameter. We have thus found a first order transition beyond the standard FFLO transition line. 
This discussion about the order of the transition is the exact analogue of the one we made in the introduction for the 
T = situation. 

Naturally it is of interest to minimize /3 since it is rather natural to expect that the states corresponding to the 
minimum will lead to the stronger instability toward the first order transition. We show now that (3 > —5, the 
equality /3 = — | being obtained for any real order parameter. We make use of : 

J [A2(VA*)2 - |A|2|VA|2] + ex. = J [AVA* - c.c.]^ < (13) 

where the equality occurs only for a real order parameter (within an irrelevant overall constant phase factor), and : 

J [A2(VA*)2 + 2|A|2|VA|2] + ex. = - J {A^AV^A* + e.c.) = 2ql J \A\^ (14) 

where the last step makes specific use of the form Eq.(ll) for the order parameter. We have then : 

- 2pql J |A|4 = J A2(VA*)2 + C.C. < J [A^{VA*f - ^[A^VA*f - |A|2|VA|2]] + e.c. (15) 

= ^ y [A\VA*f + 2|A|2|VA|2] + C.C. = ^ql J \A\' 

hence /? > — |. 

Together with /3 it is also intuitively convenient to consider 7 defined by : 

7 / dr\A\^ = q,' f dr[V\A\r = ^(qi - q^f A^,A*^^A^,A*^^ (16) 



which is easily seen to satisfy 7=1-/3 from Eq.(14). We are thus interested in maximizing 7. From the first 
expression in Eq.(16) it is intuitively clear that 7 will be large when the (properly normalized) order parameter has 
strong spatial variations. In particular it will be better for |A(r)p to have many nodes. This is more easily achieved if 
A(r) has no imaginary part since in this case one has only to require that the real part is zero. This makes intuitively 
reasonable that 7 is maximized by a real order parameter. One can also come to this conclusion from the Fourier 
expansion in Eq.(16) (or also from Eq.(12)) : it is of interest to have as often as possible opposite wavevectors so 
that (qi — €[2)^ takes as much as possible its maximum value, namely 4 . One is thus naturally led to an order 
parameter which is a combination of cos(qo.r + V'qo) with real coefficients, the directions of the qo 's being free. This 
corresponds merely to require that A(r) is real in Eq.(ll) by taking A-q^ = Aq^^. On the other hand, since any real 
order parameter gives the maximum 7, it is not necessary that these cosines have equal weight. 

We can now come back to our free energy Eq.(9) and find the best solution within the LO subspace, since we have 
found that the minimum (3 is —1/3, as soon as the order parameter is real, which implies A-q^, = Aq^ in Eq.(ll). 
We introduce a measure A of the amplitude of the order parameter by setting : 

/ dr\A\' = ^A^,Al=N2A^ (17) 

where by definition A'2 = iV is the number of plane waves coming in Eq.(ll). If all the planed waves have same 
amplitude, A is just the common value of these amplitudes. Then we define N4 and Nq by : 
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/ |A|4 = ^ Aq, A;^ Aq3A;^ = N,A^ (18) 

•' qi 

and : 

[ \Af = J2^^^KA^^KA^^Al = NeA' (19) 

For the simple plane wave considered in Eq.(8), we had N2 = N4 = Nq = 1. For the real order parameter we are 
interested in, the set of wavevectors {q^} is made of N/2 pairs. If all the planed waves have same amplitude, one finds 
| p9[ for example by simple counting that N4 — 3N{N — 1) and Nq = 5N{3N'^ — 9N + 8) . Actually, once it is recognized 
that, from the counting procedure, N4 and Nq are polynomials in N of order 2 and 3 respectively, the coefficient may 
easily be found by considering the cases N = 2,4, 6. With these notations the free energy Eq.(9) reduces to : 

F = N^A'iA, + - iV4A4[i - ^q^] + ^N^A^ (20) 

where we have set a — (for the simple plane wave considered in Eq.(8), we had a ~ 1). 

We proceed now as we have done for Eq.(8). Minimizing F with respect to we find for the extremum the 
condition 9^ = f — |a(-^4/-^2)A^, which implies that, in our considerations, we have for A^ an upper bound 
^max — {2/'ia)N2/N/i. This leads, for the value of the free energy F at this extremum, to : 

^ = ^^^(^0 - ^) + - -) + A^(^ - — ^) (21) 

Now the standard FFLO solution corresponds to — Since we are interested in a better solution we want 
v4o > which makes the first term of F/A"^ positive. On the other hand, for A = A,-aax^ the sum of the last two 
terms in the right hand side (r.h.s.) of Eq.(21) can be written as A'^^^{Ni/24a)[h{a - |)^ + 2N2Nq/NI - |]. This 
is always positive since we have N2Nq/N^ ^ 1 ( this results directly from / dr|Ap./ (ir|A|^ > [/ (ir|A|'']^). If we 
assume that the second term in the r.h.s. of Eq.(21) is positive, this implies that the free energy is always positive. 
Therefore if we want to find non positive values for F it is necessary to have a negative coefficient for the second 
term in the r.h.s. of Eq.(21), which means a < |. This is possible since the minimum /3 = — i we have found above 
corresponds to a minimum a ~ ^ 
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Then the quadratic form Eq.(21) can be zero if we meet the condition 4(Ao - ^)(2^^pa - 5a^) < (^ - 1) . This 
leads to the following result for the transition line to the FFLO phase : 

5 1 (00) 



In the (j.i — T) plane this line is higher than the standard second order FFLO transition line, which is given by Aq 
On the other hand the value of Am which gives F = at the threshold given by Eq.(21) is : 
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^4 



which is clearly nonzero and our transition is quite explicitely a first order transition (note that the condition A^ < 
Amax is necessarily satisfied since F is zero for A^ and positive for Amax )• 



Coming back to the location of the transition line Eq.(23) we consider now how to optimize it. First we see from 
the numerator that it is advantageous to lower a as much as possible. Indeed our minimum a = | is quite close to the 
limiting value a = | so the variation of the denominator with a is irrelevant. Hence we are lead to make the fourth 
order term in the free energy Eq.(9) as negative as possible by minimizing /3, as we anticipated at the beginning of 
section IV. Once we have taken a = |, we see that it is of interest to take N2Nq/N^ as small as possible (we have 
seen that it is bounded from below by 1). From their definitions Eq.(18) and Eq.(19) we can evaluate N4 and iVg for 
a general order parameter Eq.(ll) A(r) = J2 ^q; exp(jqi.r) = 2^ l^qj cos(qi.r + <Pq.). One finds : 
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7V4A4 = 3(5:iAq.Pf -Sj^lAqJ^ (24) 

i i 

and : 

iVg = 15(^ I Aq. \Y - 45(^ I Aq. J2 I Aq. |^ + 40 ^ | A^. 1^ (25) 

i i i i 

Actually one can recognize that, for symmetry reasons, the results involve only J^i l^qj" with n — 2,4,6, because 
odd powers of cosines average to zero. Hence the results assume necessarily the general form given by Eq.(24) and 
Eq.(25). Then the coefficients are easily obtained from the specific case where all the amplitudes |Aq. | are equal. 

Since the wavevectors are paired it is now more convenient to sum over pairs from now on. Defining a; = 
l^qi ?/ Y.i l^q, 1^ (implying Y,i o.i = 1) we have : 

9 6-95*2 + 4S'3 

^,N,N,lNl = ^^zV^ (26) 

where 5*2 = X^i '^i ^-nd 5*3 — X^i'^f- When we have a single pair ^2 = 6*3 = f and the r.h.s. of Eq.(26) is equal 
to 1. We show now that it is otherwise larger than f . Since (2 — S2Y < 4 — 35*2 because ^2 < 1 when we have 
more than a single pair, it is enough to prove that 3(f — 6*2) > 2(f — ^3). This is in turn verified because we can 
write for the left-hand side 1^82 — (Xi*^*)^ ~ Xi ~ 2Xi<i liOj- In a similar way we have in the right-hand 



our 



side 1-53 = = 'ill^<j a^aj{a, + aj) + QY.i<j<k (^^^^i^^k ^ SJ2^<J ataj - ^J2^<j<k (^^(^j"■k■ So 

statement is correct since it is just equivalent to Xj<j<fe '^i^-j'^k > (the equality holds when we have only two pairs 
since one can not have three different indices). Therefore we come to the conclusion that N2Nq/N2 is minimized when 
we take a single pair of plane waves with wavevectors (q, — q), corresponding to a simple order parameter proportional 
to cos(q.r), which has hence a planar symmetry. It is unfavourable to increase the number of plane waves. This is 
more easily seen in the particular situation where these N plane waves have the same amplitude. In this case we 
merely have N2Nq/N^ = [15A''^ — 45N + 40]/9(7V — f )^ which increases regularly with increasing N, and is minimum 
for TV = 2 . In this last case N2NQ/NI — It is interesting to remark that this conclusion is opposite to what one 
would obtain by considering the fourth order term alone in Eq.(9) and omitting the sixth order one. Since the fourth 
order term grows (compared to the second order one) with the number N plane waves, one would conclude that it is 
better to increase this number. The opposite turns out to be true because the sixth order term grows even faster with 
the number of plane waves. This shows quite clearly that, in contrast with what one might hope, the consideration 
of the fourth order term is not enough to conclude about the actual ground state of system. 

We find then explicitely Aq = + 2.02 fO~^. We note that this is quite close to the standard FFLO transition 
itself. For comparison the standard Clogston-Chandrasekhar |^,^ first order transition is given, as we have seen, by 
^0 = I so the difference with the standard FFLO transition is —1.39 I0~^ which is also rather small. Since we know 
that the standard Clogston-Chandrasekhar transition and standard FFLO transition stay close beyond the vicinity of 
the TCP and that this proximity between these two lines extends down to zero temperature, it seems quite possible 
that the same is true for our first order transition line. This is indeed what has been found numerically by Matsuo et 
al. . On the other hand we know from the work of Burkhardt and Rainer that the transition is second order 
in a 2D situation. So it is of some interest to consider formally an arbitrary dimension to see how one goes from the 
2D to the 3D situation. For the case of our planar order parameter, there is no difficulty in working with an arbitrary 
dimension D. As we have mentionned at the end of section II, one has just to modify some coefficients in Eq.(7) : here 
one has to multiply the terms by 3/D and the terms by 15/D{D + 2). Reproducing then the above analysis, 
we obtain that the FFLO transition occurs for Aq = -^^^ with a wavevector (f — One finds again that, in order 

to obtain a phase transition higher than FFLO, the coefficient of the second term in the r.h.s. of Eq.(2I) has to be 

negative. But this coefficient is now in general iV4( "^''^^^ — l)/4 with explicitely a = |. So we find Dc = 2.5 as the 
critical dimension to have a first order transition, which is halfway between the physical situations D = 2 and D — 3. 
Hence, in contrast to what one might think, the first order transition we have found is not 'very near to be second 
order', despite the proximity between the first order transition line and the standard second order FFLO line. The 
result Eq.(21) for the location of the transition for this planar state becomes : 

\2D 15D(7D - 10) ^ ' 

with the value of A„i of A at the threshold given by : 
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-2 4 2L>-5 

In two dimensions we can again make use of Eq.(20) to compare all the possible states. As we have seen the transition 
is second order and we can neglect the sixth order term as in the original LO analysis. The fourth order term is 
positive and, roughly speaking, below the transition the state with smaller fourth order term is selected. This leads 
to take N4/N2 = 3(2 — 5'2)/2 as small as possible. So we find again that the single pair state = 2 is selected, which 
corresponds to the state investigated by Burkhardt and Rainer . 



In closing this section, let us note that we have have found a non- degenerate minimum for the free energy, namely 
the planar state cos(q.r) corresponding to two plane waves with opposite wavevectors and equal amplitude in Eq.(ll). 
However by adding to this order parameter another cosine with a small amplitude, namely ecos(qi.r), we can have 
a free energy arbitrarily close to the planar state free energy. Conversely the location for the transition can be made 
arbitrarily close to the planar state transition. These states are highly degenerate since the relative directions of q 
and qi do not matter. These are some kind of excited states for our system with arbitrarily small energy and it is 
possible to speculate that they play a role in the physics of our system. 



V. BEYOND THE LO SUBSPACE 



In the preceding section we have restricted to the LO subspace our search for the state with lowest free energy 
at the transition and therefore with highest critical temperature. While this restriction is justified when the free 
energy expansion is limited to the fourth order term (this term being treated as a perturbation), it is no longer valid 
when the sixth order term is included. And indeed the actual minimum state does not belong to the LO subspace, 
since it can be checked that Eq.(ll) can not satisfy the Euler-Lagrange equation corresponding to Eq.(lO) because of 
the non linear terms. However we have found in the preceding section a non degenerate minimum corresponding to 
the planar state. The selection of the wavevectors of this state was due to the fourth order term while the number 
of components was controled by the sixth order one, which is otherwise quite structureless. If the nonlinearities 
produce only a rather small change on the solution, it is reasonable to look for the actual solution as being 'near' 
the solution cos(q.r) we have found in section IV. In particular it is reasonable to look for an order parameter with 
a one-dimensional dependence on r and we will also assume it to be real. And indeed we will find a solution which 
is quite near cos(q.r). Although this argumentation is only a self-consistent one (we can not exclude that a strong 
nonlinear modification of a three-dimensional solution produces the actual minimum), we note that the numerical 
exploration of Houzet et al has indeed produced a real one-dimensional order parameter for the minimum. 

With a one-dimensional real order parameter the reduced free energy simplifies into : 

; 2 1 A '2 i 1 A "21 [ \ i ^a2a'2i,1/"j^a6 



dx MoA - -A ^ + -A ^] - / da; [-A* - -A^A ^] + - / A^ (29) 
J 35J46 8J 

where A' and A" are first and second derivative of A with respect to x. Although we deal with a nonlinear problem 
we can still minimize with respect to the amplitude of the order parameter, just as we have done at the end of the 
preceding section (this works actually also for a three-dimensional order parameter). We set A(x) = a5{x), where 
5{x) is a normalized spatial function (for example by j 5^ = \ ). 5{x) gives the shape of the order parameter while 
a corresponds to its amplitude. When we substitute in Eq.(29) we obtain for F/a^ a quadratic form in . Writing 
again that, for a specific function 5{x), this form has double root at the transition leads us to an expression for Aq 
which does not depend anymore on the normalization of 5 : 

[js^-fjssr? , /K^-K^ .... 



and we want now to maximize Aq. Similarly the condition that the free energy is zero at the transition imply : 

-A2A'2 (31) 



or equivalently 
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Ao / A^ = / I A' + iA'2 - iA"2 (32) 
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which gives the ampHtude of the order parameter. 



The condition that Aq is maximum can be obtained from Eq.(30) but it is easier to deduce it from Eq.(29). One 
finds the ordinary nonhnear differential equation : 

Ia"" + (i - Ia^)A" - + ^A^ - Ia^ + AoA = (33) 

5 3 6 6 8 2 

(analyticaUy this equation can be integrated once and then reduced to a second order nonhnear differential equation but 
this is of no real help for numerics). One checks readily that the cosine form of the order parameter A(a::) = acos{qx) 
is not a solution. It satisfies only the linear part of the equation, which leads to Ao — (f /3 — Looking for the 

maximum of gives the FFLO result Aq ~ 5/36 ~ 0.1389 with Qppi^Q — 5/6. This FFLO solution corresponds to 
keep only the second term in the r.h.s. of Eq.(30). This term produces for Aq a kind of (inverted) effective potential 
g^/3 — g*/5 which has a strong maximum for the FFLO solution. When we write a Fourier expansion of the solution 
this potential selects quite effectively the wavevectors in the close vicinity of the FFLO result (Jfflo- 

Indeed the numerical exploration of Eq.(33) for Aq very close to the maximum gives for the solutions one or two 
basis frequencies qo and qi which are close to qpFLO- The other frequencies appearing in a Fourier analysis are simply 
odd combinations of go and qi like 2qo ± gi, 2qi ± go, 3go, 4gi ± goj etc. The weights of these frequencies depend on 
their order, i.e. the higher the frequency, the smaller the weight. This shows explicitely the strength of the effective 
FFLO potential in Ao (the second term in the r.h.s. of Eq.(30)). We have investigated analytically the efficiency of 
a small frequency splitting by writing qi = qo + e with e — > (but e ^ because this limit is singular). However a 
single frequency qi ~ qo turned out to be always better. We are finally lead to the conclusion that going beyond the 
LO subspace produces only small corrections due to nonlinearities. These corrections correspond to odd harmonics 
of the fundamental frequency q and they are small because of the efficiency of the effective FFLO potential. Actually 
if we assume from the start that the corrections to 6o{x) = cos{qx) are small, we can write d{x) = So{x) + Si{x), with 
di{x) small, and perform a first order expansion of Eq.(30). One finds readily that it is most favorable to take 6i{x) 
proportional to cos'^(ga;) (another contribution from cos^(ga;) is quite small), but one has to go to second order to 
find the amplitude. Actually, once it is proved that this harmonic expansion is the correct answer, it is much easier 
to minimize Eq.(30) numerically which avoids cumbersome calculations. Specifically we considered the trial function 
S{x) — cos(ga;) + a3 cos{3qx + + cos(5g + 02) and numerically maximized Ao (some more complete trial functions 
that we also tried eventually reduced to this form when maximized). We found 03 = —1.33.10"^, 05 = 1. 62.10"^, 
q = 0.793, (j>i^2 = and Aq = 0.141604. This form is very close to the cosine solution and the improvement in Ao 
compared to our result Eq.(22) for the cosine solution = 31/220 ~ 0.14091 is pretty small in absolute values, 
although it gives a significant increase of 7 10~^ to our gain of 2.02 10"'^ compared to the FFLO result. Our result 
q = 0.793 has to be compared to our result in LO subspace from Eq.(23) q = 0.829 and to qpFLO = 0.913. 



VI. CONCLUSION 



In this paper we have explored analytically the nature of the transition to the FFLO superfiuid phases in the 
vicinity of the tricritical point, where these phases begin to appear. This region is convenient for the analytical 
study we make because, in the vicinity of this point, one can make use of an expansion of the free energy up to 
sixth order, both in order parameter amplitude and in wavevector. Despite this simplification one has still a complex 
nonlinear problem to solve. We have first explored the minimization of this free energy within the LO subspace, made 
of arbitrary superpositions of plane waves. We have seen that the standard second order FFLO phase transition is 
unstable and that a first order transition occurs at higher temperature. Within this subspace we have shown that it 
is favorable to have a real order parameter. Then among these states we have shown that those with the smallest 
number of plane waves are prefered. This leads to retain only two plane waves, corresponding to an order parameter 
with a cos(qo.r) dependence, in agreement with preceding work p^ . Finally we have shown that, when releasing the 
constraint of working within the LO subspace, the order parameter at the transition is only very slightly modified by 
higher harmonics contributions and we have been able to ascribe this result to the strong selection of the wavevector 
caused by the second order terms of the free energy, corresponding physically to the standard FFLO transition. 

We are very grateful to X. Leyronas for stimulating discussions and to K. Nagai for making us aware of his work, 
during a meeting very nicely organized by G. Eska and H. Brandt for D. Rainer's 60th birthday. 
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